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ABSTRACT 

Context. Two classes of gamma-ray bursts (GRBs), short and long, have been determined without any doubts, and are usually pre¬ 
scribed to different physical scenarios. A third class, intermediate in Tgo durations has been reported in the datasets of BATSE, Swift, 
RHESSI, and possibly BeppoSAX. The latest release of > 1500 GRBs observed by Fermi gives an opportunity to further investigate 
the duration distribution. 

Aims. The aim of this paper is to investigate whether a third class is present in the log Tgo distribution, or whether it is described by a 
bimodal distribution. 

Methods. A standard fitting of a mixture of Gaussians was applied to 25 histograms with different binnings. 

Results. Different binnings give various values of the fitting parameters, as well as the shape of the fitted curve. Among five statisti¬ 
cally significant fits, none is trimodal. 

Conclusions. Locations of the Gaussian components are in agreement with previous works. However, a trimodal distribution, under¬ 
stood in the sense of having three distinct peaks, is not found for any binning. It is concluded that the duration distribution in the Fermi 
data is well described by a mixture of three log-normal distributions, but it is intrinsically bimodal, hence no third class is present in 
the Tgo data of Fermi. It is suggested that the log-normal fit may not be an adequate model. 

Key words, gamma rays: general - methods: data analysis - methods: statistical 


1. Introduction 


Gamma-ray bursts (GRBs) are most commonly classified based 
on their duration time Tgo (time during which 90% of the burst 
fiuence is accumulated, starting from the time at which 5% of 
the total fiuence was detected). [Mazets et al.| ( |1981| l first ob¬ 
served the bimodal distribution of Tgo drawn for 143 events de¬ 
tected in the KONUS experiment. |Kouveliotou et ^ ( | 1 993) 1 also 
found a bimodal structure in the log Tgo distribution in BATSE 
IB dataset; based on this data GRBs are divided into short 
(Tgo < 2s) and long (Tgo > 2 s) classes. |Horvath| ( |1998| l and 
Mukheijee et al. ( 1998| l independently discovered a third peak in 
the duration distribution in the BATSE 3B catalog ( |Meegan et al.| 
1996| l, located between the short and long peaks, and the statis¬ 
tical existence of this intermediate class was supported ( [Horvath 
2002| l with the use of BATSE 4B data. Evidence for a third com- 
ponent in log Tgo was also found in the Swift data ([Horvath et 


al.||2008a|bl Zhang & Choi 2008 Huja, Meszaros & Rfpa|20(^ 

Horvath et al. 2010 1 . Other datasets, i.e., RHESSI (Rrpa et al. 
2009 1 and BeppoSAX ( [Horvath 2009 1 , are both in agreement 
with earlier results regarding the bimodal distribution, and the 
detection of a third component was established on a lower sig¬ 
nificance level (compared to BATSE and Swift) owing to less 
populated samples. Hence, four different satellites provided hints 
about the existence of a third class of GRBs. Contrary to this, 
the durations observed by INTEGRAL have a unimodal distri¬ 
bution that extends to the shortest timescales as a power law 
([Savchenko, Neronov & Courvoisier|2012|l. 


Only one dataset (BATSE 3B) was tmly trimodal in the sense 
of having three peaks. In the others, a three-Gaussian distribution 
was found to follow the histogram better than a two-Gaussian, 
but those fits yielded only two peaks, so although statistical anal¬ 


yses support the presence of a third class, its existence is still elu¬ 
sive and may be due to selection and instmment effects, or might 
be described by a distribution that is not necessarily a mixture of 
Gaussians. The latest release is due to Fermi GBM observations 
( |Grub^ et al.[2014][von Kienlin et al.[2014[ ) and consists of 1568 
GRB^with calculated Tgo. This sample, to the best of the au¬ 
thor’s knowledge, has not yet been investigated for the presence 
of an intermediate class; only [Horvath et aL] ( [2012[ l and [Qin ^ 
al. (2013|l conducted research on a subsample consisting of 425 


bursts from the first release of the catalog. 

Previous analyses, although consistent with each other (ex¬ 
cept INTEGRAL), showed variances in locations of the Gaus¬ 
sian components, their dispersion, relative heights, and statisti¬ 
cal significance of a trimodal fit compared to a bimodal distri¬ 
bution. The aim of this article is to perform a statistical anal¬ 
ysis of the Fermi database in order to verify the existence of 
the intermediate-duration GRB class. The current sample con¬ 
stitutes ~ 83% of long GRBs (Tgo > 2 s), which is a larger per¬ 
centage than in previous catalogs. A conventional fitting is 
applied. Since BATSE 4B (Horvath 20021, the maximum log- 
likelihood method has been used. In the case of the Swift data 
this was fully justified because of the small size of the popula¬ 
tion ( Horvath et al.[200^ b I. However, in the first analysis of the 
BATSE 3B data ( [Horvath 1998| l, which consisted of 797 GRBs, 
three classes were detected with the;^^^ method, which was later 
supported by the maximum log-likelihood method on a bigger 
dataset ( [Horvath|2002[ l. Hence, the number of GRBs under con¬ 
sideration in this study should not be too small to conduct a con¬ 
ventional fitting analysis. To verify whether the fitting is ap- 


* http://heasarc.gsfc.nasa.gov/W3Browse/fermi/ 
fermigbrst.html accessed on March 12, 2015. 
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plicable, the BATSE 4B data (2041 GRBs) were re-examined 
and results similar to those of |Horvafli| ( |2002| l were obtained. To 
be sure that the method chosen does not rule out any feature in 
the data under consideration, the maximum log-likelihood was 
also applied, but since the results were not different from those 
obtained using the fitting, they are not reported here or men¬ 
tioned in what follows. 

This article is organized as follows. Section 2 presents the;^^^ 
fittings of different models. In Sect. 3 a comparison with previ¬ 
ous results is discussed. Section 4 is devoted to discussion, and 
Sect. 5 gives concluding remarks. The computer algebra system 
Mathematica® vlO.0.2 is applied throughout this paper. 


2. Jhex^ fittings 

2.1. Sample 

In this section a dataset of 1566 GRBs is used. Two durations 
(the shortest and the longest) were excluded because, for every 
binning, they were separated from the rest of the distribution by 
empty bins. 


2.2. Standard normal distributions 

A least squares fitting of a multi-component log-normal distribu¬ 
tion to the dataset of durations Tgo is performed, i.e., a mixture 
of Gaussians 


k 


k 


/=1 /=1 


A, 

-exp 

V^CT; 


2cr] , 


( 1 ) 


is fitted to a histogram of log Tgo. A significance level of 
a - 0.05 is adopted. Twenty-five binnings are applied, defined 
by the bin widths w from 0.30 to 0.06 with a step of 0.01. The 
corresponding number of bins ranges from 15 to 69. A number 
of binnings is chosen rather than a neutral bin width established 
according to some conventional rule (e.g., Freedman-Diaconis, 
Scott, Knuth) because of |Koen & Bere et al.| ( |2012l l, who found 
two fits to Swift data that could not be rejected, hence restricting 
the analysis to only one binning might be concealing. For each 
binning the same fitting procedure is performed as described. 

The of a fit is calculated in a standard way as 


A 


2 


Z 


{Oi - Elf 
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where G, is the observed /-th value, is the value expected 
based on the fit, and N is the number of bins. The number of de¬ 
grees of freedom, dof, of the;^f^ test statistic is dof - N - m- 1, 
where m is the number of parameters used; for a k-Gaussian, 
m - 3k. First, a single-Gaussian fit is performed for all binnings. 
The;^f^ values range from 12 180 to 54548, with p-values being 
numerically equal to zero in all cases. It follows from the huge 
X^ values that it is extremely unlikely that the single-Gaussian fit 
describes the log Tgo distribution correctly. 

Next, a two-Gaussian fit is performed and the resulting x^ 
are much lower; however, the majority of p-values indicates that 
the distributions do not follow the data well at the a = 0.05 sig¬ 
nificance level. Fitted parameters are gathered in Table with 
p-values greater than a written in bold, while the curves are dis¬ 
played in Fig. The bimodal structure is well represented by a 
two-Gaussian fit in four out of five cases displayed in Fig.[2 The 
statistically unsignificant fit for w = 0.25 is also shown due to a 
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Fig. 1. Two-Gaussian fits (solid black) to log Tgo distributions. All ex¬ 
cept w = 0.25 are statistically significant. Colored dashed curves are the 
components of the mixture distribution. 
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Table 1. Parameters of a two-Gaussian fit 


w 

i 

Fi 

o-i 

A,- 


p-val 

0.27 

1 

-0.042 

0.595 

100.2 

6.467 

0.692 

2 

1.477 

0.465 

325.7 

0.26 

1 

-0.063 

0.569 

92.30 

12.23 

0.270 

2 

1.475 

0.473 

318.8 

0.25 

1 

-0.125 

0.510 

79.55 

22.00 

0.024 

2 

1.453 

0.494 

316.1 

0.20 

1 

-0.049 

0.611 

73.34 

22.09 

0.077 

2 

1.473 

0.476 

243.6 

0.13 

1 

-0.030 

0.607 

48.56 

33.74 

0.142 

2 

1.480 

0.468 

157.2 


Table 2. Parameters of a three-Gaussian fit 


W 

i 

Fi 

o-i 

A,- 


p-val 


1 

-0.030 

0.603 

102.0 



0.27 

2 

1.466 

0.455 

317.1 

5.333 

0.502 


3 

2.027 

0.201 

6.014 




1 

-0.210 

0.461 

71.48 



0.26 

2 

1.119 

0.450 

128.6 

6.819 

0.448 


3 

1.598 

0.421 

208.4 




1 

-0.137 

0.492 

77.73 



0.25 

2 

1.414 

0.480 

300.4 

13.52 

0.095 


3 

1.939 

0.128 

14.49 




1 

-0.204 

0.493 

57.30 



0.20 

2 

1.221 

0.488 

144.0 

17.71 

0.087 


3 

1.665 

0.396 

113.1 




1 

-0.058 

0.581 

46.61 



0.13 

2 

1.453 

0.464 

153.7 

29.42 

0.167 


3 

1.903 

0.092 

4.328 




comparison with a three-Gaussian fit to be performed further on. 

Finally, a three-Gaussian fit is performed in the same man¬ 
ner. The resulting are naturally lower than they were in the 
previous case. Again, a majority of the corresponding p-values 
that are smaller than a indicates that in general the duration dis¬ 
tribution is not well described by a mixture of three log-normal 
components. Parameters of the five fits with p-values > a are 
gathered in Table and the fitted curves are displayed in Fig.[^ 
It is important to note that these five fits are all statistically sig¬ 
nificant, including the fit for w - 0.25, for which a hypothesis 
that the histogram is well described by a two-Gaussian fit was 
rejected. 

The probability that the third component originates from 
statistical fluctuations may be estimated by comparing a three- 
Gaussian with a two-Gaussian fit, based on an observation ([Band] 
let al.|1997| l that 

Ax^^xI-xI-xHm, (3) 

where Av = V 2 - vi is the difference in the degrees of freedom of 
fits under consideration, equal to three when Ak = 1 (compare 
with Eq. ([^), and = denotes equality in distribution. 

To elaborate to what degree a three-Gaussian fit is bet¬ 
ter than a two-Gaussian, Eq. Q is applied and the p-values 
are inferred from a x^ distribution with three degrees of free¬ 
dom. A small p-value indicates a small probability that the two- 
Gaussian model alone describes the data in comparison with a 
three-Gaussian. The results, gathered in Tableindicate that in 
all binnings there is a > 20% probability (exceeding 50% for 



logTgo 

Fig. 2. The same as Fig[^ but for a three-Gaussian fit. 
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Table 3. Improvements of a three-Gaussian over a two-Gaussian fit 


w 

Ax' 

p-value 

0.27 

1.134 

0.767 

0.26 

5.411 

0.144 

0.15 

8.480 

0.037 

0.20 

4.380 

0.223 

0.13 

4.320 

0.229 


the remaining three statistically significant cases) that a three- 
Gaussian is a significant improvent over a two-Gaussian fit. This 
does not mean that the three-Gaussian models are a good de¬ 
scription of the data, as the from Table lead to rejection 
of the null hypothesis for only one out of five binnings with p- 
values greater than a in the three-Gaussian case. For this un¬ 
significant fit, the p-value computed from means that a bad 
three-Gaussian is better than a bad two-Gaussian. For the fits that 
are statistically significant (when fitting both a two- and a three- 
Gaussian) the conslusion is that a null hypothesis cannot be re¬ 
jected for both, but a three-Gaussian describes the data better 
with a probability > 50%. Nevertheless, this probability is insuf¬ 
ficient to claim a detection; in ( |Horvath|2009| ) it was concluded 
that even a 96% significance level is too small to be considered 
evidence. 

The Fermi distribution is dominated by long GRBs, which 
constitute ~ 83% of the sample, that is manifested through 
a significantly higher dispersion of the short GRB distribu¬ 
tion. Moreover, it was proposed that the distribution for short 
GRBs was nearly flat for Tgo < 2s ( [Bromberg, Nakar & Pi- 
201 1[ Bromberg et al. 2013|l (compare with (Savchenko, 


ran 


Neronov & Courvoisier|2012|l); however, for smaller bin widths 


a statistical noise starts to dominate. This is supported by a 
very small p-value of ~ 10“^ - 10“^ for the smaller bin widths, 
w = 0.11 - 0.06, which indicates that at these binnings statistical 
noise is dominating. 

Three out of five statistically significant fittings 
(w = 0.27, 0.25, 0.13) located the third component beyond 
the main peak for long GRBs. One fit (w = 0.26) showed an 
excess at Tgo a; 10 s that might be assigned to an intermediate 
class; however, the fitted curve is very similar to a two-Gaussian. 
The last fit, with w - 0.20, detected two components (in addition 
to a peak related to short GRBs) of approximately the same 
height and comparable standard deviations in a region of long 
GRBs. For these last two fits, the dispersion of a corresponding 
peak of a two-Gaussian is greater than the dispersion of both the 
components in the three-Gaussian, hence it gives a hint toward 
the bimodality of the GRBs ( Schilling, Watkins & Watkins | 
2002|l. This is in agreement with a very high probability that 


As short GRBs were underpopulated, the three-Gaussian fit, de¬ 
spite following the observations better than a two-Gaussian, con¬ 
sisted of a component with a high dispersion. Overall, classifi¬ 
cation analyses found three components in all four samples. In¬ 
terestingly, the dataset from INTEGRAL ([Savchenko, Neronov 


& Courvoisier 


20121 yields a unimodal distribution. The latest 


the third group in a three-component model is a statistical 
fluctuation. 


3. Comparison with former resuits 

A detailed comparis on between BATSE 4B and Sw ift catalogs 
has been conducted ( Huja, Meszaros & Rfpa|[2009 1 and the re¬ 
sults were found to be consistent, also in the means of differences 
between the instru ments. The RHES SI observations were also 
taken into account (Rfpa et al. 2009 1 and roughly the same dis¬ 
tribution as in previous catalogs was reported. The BeppoSAX 
data were of a relatively low population and the analysis showed 
the presence of three components on a significance level lower 
than in previous catalogs, but two classes - the intermediate 
and long - were detected with high certainty ([Horvath[[2009|l. 


well-populated sample, based on the Fermi data, was investi¬ 
gated some time ago ( Horvath et al.|2012 l using PC A and mul¬ 
ticlustering analysis, and a three-group structure was found in a 
multidimensional parameter space including duration, total flu- 
ence, hardness ratio, and peakflux256. Hence, data from five 
satellites supported a three-component distribution of GRBs by 
means of statistical significance. 

In Eig. all of the components’ locations found by the 
above-mentioned univariate analyses are plotted. Results from 
this work, i.e., locations of the three-Gaussian components found 
by a standard log-normal fitting, are consistent with previous re¬ 
sults by means of locations as well as relative separations - mean 
durations of the five statistically significant fits are centered at 
0.745 s, 21.61s, and 67.05 s for short, intermediate, and long 
GRBs, respectively. The current sample of 1566 GRBs com¬ 
prises 17% of short GRBs and 83% of long ones, based on the 
conventional classification S 2 s. Because of an insufficient sep¬ 
aration of the intermediate and long components, it is impossi¬ 
ble to conclude what the population of intermediate class GRBs 
in Fermi data is, as the distribution is apparently bimodal and 
shows no evidence for the third class being present, and one can 
associate the intermediate and long groups with a single class. 
It is important to note that in Fermi the sensitivity at very soft 
and very hard GRBs was higher than in BATSE ( [Meegan et al.[ 
[2009| l. Soft GRBs are intermediate in duration, and hard GRBs 
have short durations. Hence, an increase in intermediate GRBs 
relative to long ones might be expected as a consequence of im¬ 
proving instruments, yet the third class remains elusive. Swift is 
more sensitive in soft bands than BATSE was, hence its dataset 
has a low fraction of short GRBs. Therefore, the group popula¬ 
tions inferred from Fermi observations are reasonable consider¬ 
ing the characteristics of the instruments. 

4. Discussion 

Among the 25 three-Gaussian fittings performed, five 
(w = 0.27, 0.26, 0.25, 0.20, 0.13) turned out to be statisti¬ 
cally significant, with p-values exceeding the significance level 
a - 0.05. Locations of the respective groups for the five fits 
are close to each other with means pi - -0.128 + 0.082 (short 
GRBs), p 2 = 1-335 + 0.156 and p^. = 1-826 + 0.185 (long 
GRBs), the error being the standard deviation of the average 
(compare with Eig.[^. Previous works on datasets from BATSE 
(|Horvath[[1998[ |2002|l an d Swift ( [Horvath et al.[[^08a[ [Huja,[ 
Meszaros & Rfpa 2009 1 indicated that a three-Gaussian is a 
better fit than a corresponding two-G aussian. On the other 
hand, a three-Gaussian fit to RHESSI (fopa et al. 20091 data 
yielded only a 93% probability of being correct compared 
to a two-Gaussian, meaning that there is a remarkable 7% 
probability that the log Tgo is well described by a two-Gaussian, 
while for BeppoSAX ( [Horvath|[2009[ ) the goodness-of-fit was 
not reported (only the maximum log-likelihoods). 

Moreover, the two greater means of the three-Gaussian from 
Eig. [^ (w = 0.20) do not satisfy t he criterion of bimodality 
( [Schilling, Watkins & Watkins[2002) l 


\llA - IIb\ > (o-A +0 -b)-S 


CT- 




.2 \ 


(4) 
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Fig. 3. Locations of the short (green squares), intermediate (blue 
crosses) and long (red triangles) GRBs from previous research and this 
work. 


with (Ta < (Tb, and S (r) being the separation factor equal to 
0.98 in the case of w = 0.20 (assuming equal weights; arbitrary 
mixtures yield an even higher value of the factor) for cr 2 and crj, 
from Table hence long GRBs (Tgo > 2 s) do not have a bi- 
modal distribution and are described by a single peak, meaning 
that the distribution over the whole range of available durations 
Tgo is bimodal, with peaks corresponding to short and long GRB 
populations. The fit for w - 0.26 does not fulfil this condition ei¬ 
ther, where the appropriate S{r) - 1.16 is taken from ([Schilling 


Watkins & Watkins |2002|l. The remaining three cases, although 


the shoulder is prominent, are also apparently bimodal. 

The relative improvements (Table indicate an enourmous 
probability, ranging from 14% to 77%, that the third component 
in a three-Gaussian fit is a chance occurrence compared to a two- 
Gaussian. 

Finally, among previous research, only the BATSE 3B data 
revealed a truly trimodal logTgo distribution, i.e., having three 
local maxima, which was not present in the following release, 
BATSE 4B (where only a bump was present), and is non-existent 
in the current BATSE catalog. Swift also observed two local max¬ 
ima, although with a prominent shoulder on the left side of the 
long GRB peak ([Horvath et al.|[2008a[), detected by means of 


the maximum log-likelihood method. Huja, Meszaros & Rfpa 


( [2009| l also obtained a bimodal distribution with a bump on one 
side of the long GRB peak, although somewhat weaker. The ex¬ 
planation may be that for the sample of 388 GRBs, a maximum 
log-likelihood method is more robust than applying a fitting. 
The latter may be a drawback in undersampled populations. Eor- 
tunately, this is not a case i n the Fermi data. The RHESSI dis¬ 


tribution (Rfpa et al. 


fit with a shoulder, w 


20091 is also characterized by a bimodal 


lile in BeppoSAX clearly separated inter¬ 


mediate and long groups were found (|Horvath|2009|l; however, 


no short GRBs were detected in the log Ego distribution there. 
This is most likely due to a low trigger efficiency to short GRBs, 
which are highly underpopulated in that sample. Interestingly, an 
early analysis of 222 GRBs from Swift ( [Horvath et al.|2008b| de¬ 
tected the short class, while the long GRBs were unimodal, again 
with a bump. Einally, the standard three-Gaussian fits passed the 
Anderson-Darling test, hence a mixture of log-normal distribu¬ 
tions account well for the observed durations. 

The distributions fitted in this paper are strongly dependent 
on the binning applied when the locations of components and 
their relative amplitudes are considered. No statistically signifi¬ 
cant trimodal fit was found, although a shoulder (in three out of 
five fits) was detected beyond the region where previous works 
found the intermediate class, which is a surprising result. Still, 
a three-Gaussian fit is a better fit than a two-Gaussian, statisti¬ 
cally speaking. However, it is arguable whether this confirms the 
existence of the third class of GRBs. As the sum of two normal 
distributions is skewed, which is the apparent case in the Fermi 
sample, the underlying multimodal distribution may not neces¬ 
sarily be composed of Gaussians. Recently, [Zitouni et al.| ( [2015] l 
suggested that the duration distribution corresponding to the col- 
lapsar scenario (associated with long GRBs) might not necessar¬ 
ily be symmetrical because of a non-symmetrical distribution of 
envelope masses of the progenitors. 

Moreover, the more components the mixture distribution 
has, the more parameters are available for the curve to fit the 
histogram. This may be the primary reason for the calculated 
goodness-of-fit. On the other hand, in terms of statistical ac¬ 
curacy, the hypothesis that three fundamental components are 
needed to describe observed durations is corroborated. The true 
underlying form of the distribution remains obscure. Short GRBs 
have been underpopulated in observations ever since, and some 
instrument biases have been proposed to account for the emer¬ 
gence of a third peak. The trimodality in the BATSE 3B cat¬ 
alog was greatly diminished in the 4B version, which may be 
directly attributed to collecting a more complete sample. The 
Fermi database contains 75% of the number of GRBs observed 
in BATSE current catalog, hence further observations, as well as 
theoretical models that could account for the diversity of GRB 
events in more detail may clarify the properties of the Ego distri¬ 
bution and determine whether an intermediate class of GRBs is 
a physical or a statistical phenomenom. 

Einally, a mixture of log-normal Gaussians may not be a 
proper model for the Ego distribution, and a mixture of intrin¬ 
sically skewed distributions may be a better explanation of the 
observed features of the histogram. On the other hand, the dura¬ 
tion distribution might not be a sum of components defined on 
(-CX), -i-oo), but might be a piecewise function. 


5. Conclusions 

1. A mixture of three standard normal distributions was 
found to be statistically significant in describing the 
log Ego distribution for 1566 Fermi GRBs for bin widths 
w = 0.27, 0.26, 0.25, 0.20, 0.13. Average locations of the 
components are equal to 0.745 s, 21.61s, and 67.05 s for 
short, intermediate, and long GRBs, respectively, or -0.128, 
1.335, and 1.826 in log-scale. These results are in agreement 
with values obtained from previous catalogs: BATSE, Swift, 
RHESSI, and BeppoSAX. 

2. The relative improvements of a three-Gaussian fit over a two- 
Gaussian imply that there is a significant probability, varying 
from 14% to 77% among the fits, that the third component 
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is a chance occurrence. Therefore, the third GRB class is un¬ 
likely to be present in the Fermi data. 

3. None of the fits is trimodal (in the sense of having three dis¬ 
tinct peaks). Although three out of five fits show a prominent 
shoulder on the right side of the long GRB peak, the evidence 
is not sufficient to claim detection of a third class in the Fermi 
data. Therefore, the distribution of Fermi durations is intrin¬ 
sically bimodal, hence no evidence for an intermediate class 
of GRBs has been found. 

4. The observed asymmetry may come from an underlying dis¬ 
tribution composed of two skewed components, or it may be 
a piecewise function. 
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